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'. Abstract 

We analyze (3+l)-dimensional black-hole space-times in spontaneously broken 
CN . Yang-Mills gauge theories that have been recently presented as candidates for an eva- 

sion of the scalar-no-hair theorem. Although we show that in principle the conditions 
O ! ioT the no-hair theorem do not apply to this case, however we prove that the 'spirit' 

^ ' of the theorem is not violated, in the sense that there exist instabilities, in both the 

sphaleron and gravitational sectors. The instability analysis of the sphaleron sector, 



O 



•n , which was expected to be unstable for topological reasons, is performed by means of 

^ ' a variational method. As shown, there exist modes in this sector that are unstable 

against linear perturbations. Instabilities exist also in the gravitational sector. A 
Q_i' method for counting the gravitational unstable modes, which utilizes a catastrophe- 

^ ■ theoretic approach is presented. The role of the catastrophe functional is played by 

the mass functional of the black hole. The Higgs vacuum expectation value (v.e.v.) 
is used as a control parameter, having a critical value beyond which instabilities 
^ I are turned on. The (stable) Schwarzschild solution is then understood from this 

d ' point of view. The catastrophe-theory appproach facilitates enormously a universal 

stabihty study of non-Abelian black holes, which goes beyond linearized perturba- 
tions. Some elementary entropy considerations are also presented that support the 
catastrophe theory analysis, in the sense that 'high-entropy' branches of solutions 
are shown to be relatively more stable than 'low-entropy' ones. As a partial result of 
this entropy analysis, it is also shown that there exist logarithmic divergencies in the 
entropy of matter (scalar) fields near the horizon, which are up and above the linear 
divergencies, and, unlike them, they cannot be absorbed in a renormalization of the 
gravitational coupling constant of the theory. The associated part of the entropy 
violates the classical Bekenstein-Hawking formula which is a proportionality relation 
between black-hole entropy and horizon area. Such logarithmic divergencies, which 
are associated with the presence of non-abelian gauge and Higgs fields, persist in 
the 'extreme case', where linear divergencies disappear. 
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1 Introduction 



The surprising discovery of the Bartnik-McKinnon (BM) non-trivial particle-hke 
structure ^ in the Einstein- Yang-Mills system opened many possibilities for the 
existence of non-trivial solutions to Einstein-non-Abelian-gauge systems. Indeed, 
soon after its discovery, many other self-gravitating structures with non-Abelian 
gauge fields have been discovered . These include black holes with non-trivial hair, 
thereby leading to the possibility of evading the no-hair conjecture The physical 
reason for the existence of these classical solutions is the 'balance' between the non- 
Abelian gauge-field repulsion and the gravitational attraction. Such a balance allows 
for dressing black hole solutions by non-trivial configurations (outside the horizon) 
of fields that are not associated with a Gauss-law, thereby leading to an 'apparent' 
evasion of the no-hair conjecture. 

Among such black-hole solutions, a physically interesting case is that of a sponta- 
neously broken Yang-Mills theory in a non-trivial black- hole space-time (EYMH) 0]. 
This system has been recently examined from a stability point of view, and found 
to possess an instability , thereby making the physical importance of the solution 
rather marginal, but also indicating another dimension of the no-hair conjecture, 
not considered in the original analysis, that of stability. 

In this article, we shall give more details of this stability considerations by extend- 
ing the analysis to incorporate counting of the unstable modes, and going beyond 
the linear-stability case by employing catastrophe theory Q in order to analyse 
instabilities in the gravitational sector of the solution. Catastrophe theory is a pow- 
erful mathematical tool to study or explain a variety of change of states in nature, 
and in particular a discontinuous change of states that occurs eventually despite a 
gradual (smooth) change of certain parameters of the system. In the case at hand, 
the catastrophe functional, which exhibits a discontinuous change in its behaviour, 
will be the mass of the black hole space time, whilst the control parameter, whose 
smooth change turns on the catastrophe at a given critical value, will be the vacuum 
expectation value (v.e.v.) of the Higgs field. The advantage of using the v.e.v. of 
the Higgs field as the control parameter, rather than the horizon radius as was done 
in , is that it will allow us to relate the stability of the EYMH black holes to that 
of the Schwarzschild solution, which is well known to be stable. The type of catas- 
trophe encountered will be that of a fold catastrophe. The catastrophe-theoretic 
approach allows for a universal stability study of non-abelian black hole solutions 
that goes beyond linearised perturbations; the particular use of the Higgs v.e.v. as 
a control parameter in the case of the EYMH systems allows an exact counting of 
the unstable modes. 

As part of our analysis, we shall make an attempt to associate the above catastrophe- 
theoretic considerations with some 'thermodynamic/information-theoretic' aspects 
of black hole physics, and in particular with the entropy of the black hole. By 



computing explicitly the entropy of quantum fluctuations of (scalar) matter fields 
near the horizon we shall show that 'high-entropy' branches of the solution possess 
less unstable modes (in the gravitational sector) than the 'low-entropy' ones. As a 
partial, but not less important, result of this part of our analysis, we shall also show 
that the entropy of the black hole possess linear and logarithmic divergencies. The 
linear divergencies do not violate the Bekenstein-Hawking formula relating entropy 
to the classical horizon area. The only difference is the divergent proportionality 
factors in front, which, however, can be absorbed in a conjectured renormalization 
of Newton's constant in the model This is not the case with the logarithmic 

divergencies though. The latter persist even in 'extreme black-hole' cases, where the 
linear divergencies disappear. They clearly violate the Bekenstein-Hawking formula. 
In our case they owe their presence to the non-Abelian gauge and Higgs fields. The 
presence of logarithmic divergencies in black hole physics has been noted in ref . , 
but only in examples involving truncation from (3+l)-dimensional space-times to 
(1-1-1) dimensions, and in that reference their presence had been attributed to 
this bad truncation of the four-dimensional black hole spectrum. Later on, how- 
ever, such logarithmic divergencies have been confirmed to exist in string-inspired 
dilatonic black holes in (3-1-1) dimensions 0. Their presence in our EYMH system, 
and in general in non-Abelian black holes as we shall show, indicates that such log- 
arithmic divergencies are generic in black hole space-times with non- conventional 
hair, and probably indicates information loss, even in extreme cases, associated with 
the presence of space time boundaries. This probably imphes that the entropy of 
the black hole is not only associated with classical geometric factors, but is a much 
more complicated phenomenon related to information carried by the various (inter- 
nal) black hole states. The latter phenomenon could be associated with, and may 
be offer ways out of, the usual difficulties of reconciling quantum mechanics with 
canonical quantum gravity. 



The structure of the article is as follows. In section 2 we shall discuss the no hair 
conjecture for black holes space-times with non-trivial scalar field configurations by 
following a modern approach due to Bekenstein [T^. We shall show that the proof of 
the no-hair theorem fails for the case of the EYMH system, in accordance with the 
explicit solution found in ref. IQ. In section 3 we shall present a stability analysis 
of the system based on linear perturbations. We shall demonstrate the existence of 
instabilities in the sphaleron sector, following a variational approach which is an ex- 
tension of the approach of Volkov and Gal'tsov |n[] to study particle-like solutions. 
We shall also present arguments for counting the unstable modes in the sphaleron 
sector of the theory. In section 4 we shall present a method for counting the unsta- 
ble modes in the gravitational sector by going beyond the linearised-perturbation 
analysis using catastrophe theory, with the mass functional of the black hole as the 
catastrophe functional and the Higgs v.e.v. as the control parameter. In section 5, 
in connection with the latter approach, we shall estimate the entropy of the vari- 
ous branches of the solution using a WKB approximation. We shall show that the 
high-entropy branch of solutions is relatively more stable (in the sense of possesssing 



fewer unstable modes) than the low-entropy branch. As we have already mentioned, 
we shall also discuss the existence of logarithmic divergencies in the entropy, associ- 
ated with the presence of non-trivial hair in the black hole, in certain extreme cases, 
and we shall argue about an explicit violation of the classical Bekenstein-Hawking 
entropy formula, indicating a different (information oriented) role of the black hole 
entropy. Conclusions and outlook will be presented in section 6. Some technical 
aspects of our approach will be discussed in two appendices. 



: Bypassing Bekenstein's no-hair theorem in 
EYMH systems 

Recently, Bekenstein presented a modern elegant proof of the no-hair theorem for 



black holes, which covers a variety of cases with scalar fields ||T0[- The theorem is 
formulated in such a way so as to rule out a multicomponent scalar field dressing an 
asymptotically flat, static, spherically-symmetric black hole. The basic assumption 
of the theorem is that the scalar field is minimally coupled to gravity and bears a non- 
negative energy density as seen by any observer, and the proof relies on very general 
principles, such as energy-momentum conservation and the Einstein equations. From 
the positivity assumption and the conservation equations for the energy momentum 
tensor Tmn of the theory, V^'^Tmn = 0, one obtains for a spherically-symmetric 
space-time background the condition that near the horizon the radial component of 
the energy-momentum tensor and its first derivative are negative 

t; < 0, (t;)' < (i) 

with the prime denoting differentiation with respect to r. This implies that in such 
systems there must be regions in space, outside the horizon where both quantities 
in (p change sign. This contradicts the results following from Einstein's equations 
though [0, and this contradiction constitutes the proof of the no-hair theorem, 
since the only allowed non-trivial configurations are Schwarzschild black holes. We 
note, in passing, that there are known exceptions to the original version of the no- 
hair theorem , such as conformal scalar fields coupled to gravity, which come from 
the fact that in such theories the scalar fields diverge at the horizon of the black 



hole 12 . 



The interest for our case is that the theorem rules out the existence of non-trivial 
hair due to a Higgs field with a double (or multiple) well potential, as is the case for 
spontaneous symmetry breaking. Given that stability issues are not involved in the 
proof, it is of interest to reconcile the results of the theorem with the situation in our 
case of EYMH systems, where at least we know that an explicit solution with non- 
trivial hair exists albeit unstable 0. As we shall show below, the formal reason 
for bypassing the modern version of the no-hair theorem [|l^ lies in the violation 
of the key relation among the components of the stress tensor, T/ = Tg, shown to 



hold in the case of ref. |jl0[. The physical reason for the 'elusion' of the above no- 
hair conjecture lies in the fact that the presence of the repulsive non-Abelian gauge 
interactions balance the gravitational attraction, by producing terms that make the 
result (0) inapplicable in the present case. Below we shall demonstrate this in a 
mathematically rigorous way. 

To this end, consider the EYMH theory with Lagrangian 

Ceymh = \^^\Fmn? + \(t>MM? + \\dM<P? + V{(t>)] (2) 

where Am denotes the Yang-Mills field, Fmn its field strength, is the Higgs field 
and V{(t)) its potential. All the indices are contracted with the help of the background 
gravitational tensor Qmn- In the spirit of Bekenstein's modern version of the no-hair 
theorem, we now examine the energy- momentum tensor of the model (0). It can be 
written in the form 

SttTmJV = -£9MN + ^ [fmpFn'' + '^AmAn + ^Af^^TV^j (3) 



with £ = —C 



EYMH- 



Consider, now, an observer moving with a four-velocity u . The observer sees a 
local energy density 

P = ^ + ^ I^u^'FmpF^V + ^{u^'AmY + {u^'dMC^r"^ , u'^UM = -1. (4) 

To simplify the situation let us consider a space-time with a time-like killing vec- 
tor, and suppose that the observer moves along this killing vector. Then u^'^dM<P = 
and by an appropriate gauge choice u^^Am = = u^'^Fmn- This gauge choice is 
compatible with the spherically-symmetric ansatz for Am of ref. [Q. Under these 
assumptions, 

p = £ (5) 

and the requirement that the local energy density as seen by any observer is non- 
negative implies 

£ > 0. (6) 



We are now in position to proceed with the announced proof of the bypassing of 
the no-hair theorem of ref. for the EYMH black hole of ref. @]. To this end 
we consider a spherically-symmetric ansatz for the space-time metric Qmn, with an 
invariant line element of the form 



r(r), A = A(r). (7) 



To make the connection with the black hole case we further assume that the space- 
time is asymptotically flat. 



From the conservation of the energy-momentum tensor, following from the invari- 
ance of the effective action under general co-ordinate transformations, one has for 
the r-component of the conservation equation 



IT 

or 



MN 







(8) 



with the prime denoting differentiation with respect to r. The spherical symmetry 
of the space time implies that Tq = T^. Hence, (j^) can be written as 



(e 2 r'^T^y - -e 2 



r'T/ + a't; + -T^ 



0. 



(9) 



Observing that the terms containing A cancel, and integrating over the radial coor- 
dinate r from the horizon to a generic distance r, one obtains 



— _ / dre'^r 
2r^ Jru 



T'Tl + -T| 
r 



(10) 



Note that the assumption that scalar invariants, such as TmnT^^ are finite on the 
horizon (in order that the latter is regular), implies that the boundary terms on the 
horizon vanish in (ITC 



It is then straightforward to obtain 



ft; + -Tq 

r 



_ r 

e 2 r 2\trrr 



Next, we consider Yang-Mills fields of the form 

A= [l + uj{r))[-T^de + feis.in6dif] 



(12) 



where r,, i = r,9,ip are the generators of the SU{2) group in spherical-polar coor- 
dinates. Ansatz ([T^) yields 



rpr 

r 

rpO 

J- a 



with (see Appendix A for details of the relevant quantities). 



-A 
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Air 



2uj 



12 



+ 



+ 



(13) 



(14) 



Substituting in (|T3|) yields 



e"2 r 



r2 



(Tn'ir) 



where S is expressed as 



An 



t\2 



+ 



2n2 



2r4 



+ 



ear 2 ^ 



b^il + u? 1 



A 



(15) 
(16) 

(17) 



We now turn to the Einstein equations for the first time, following the analysis of ref. 
[|T^. Our aim is to examine whether there is a contradiction with the requirement 
of the non-negative energy density. These equations read for our system 



SttT;. 



Integrating out the first of these yields 



r Jm r 



(19) 



where M., is a constant of integration. 



The requirement for asymptotic flatness of space-time implies the following asymp- 
totic behaviour for the energy-density functional £ ~ 0{r~^) as r cxo, so that 
A ~ 0{r~^). In order that — oo at the horizon, r — > r/j, M., is fixed by 



M, 



The second of the equations (0) can be rewritten in the form 



— — — 2 / 2 — \/ 



4v]-rT'' H 



A 3 
2r 



(20) 



(21) 



Consider, first, the behaviour of TJ^ as r — > oo. Asymptotically, 62 — >• 1, and so 
the leading behaviour of {TJ^)' is 

2, 



(22) 



We, now, note that the fields lu and have masses | and /i = ^/Xv respectively. 
From the field equations and the requirement of finite energy density their behaviour 
at infinity must then be 



uj(r) ~ — 1 + ce 2' 



(r) 



V + ae 



(23) 



for some constants c and a. Hence, the leading asymptotic behaviour of J and is 



1 

47r 



4^2 
2^ 



4r2 



(24) 



since e — > 1 asymptotically. 



The leading behaviour of (T^)', therefore, is 



'T\l 



1 

47r 



2c2 



(25) 



There are two possible cases: (i) 2v^yU > v (corresponding to A > 1/8); in this case 
i^D' > asymptotically, (ii) 2\/2jj, < v (corresponding to A < 1/8) ; then, (T^)' < 
asymptotically. 

Since JT" vanishes exponentially at infinity, and S ~ 0[r^'^] as r ^ oo, the integral 
defining T^{r) converges as r ^ oo and \T^\ decreases as r~^. 

Thus, in case (i) above, is negative and increasing as r — >• oo, and in case (ii) 
is positive and decreasing. 

Now turn to the behaviour of TJ^ at the horizon. When r c:^ Vh , S and JT" are both 
finite, and T' diverges as l/(r — r/j). Thus the main contribution to as r ~ r/^ is 



Jr. 2 



(26) 



which is finite. 

At the horizon, = 0; outside the horizon, e'" > . Hence F' > sufficiently 
close to the horizon, and, since £^ > 0, < for r sufficiently close to the horizon. 

Since JF ~ 0[r — r^] at r ^ r^, (T^)' is finite at the horizon and the leading 
contribution is 

m'(r.)^-^^+-X (27) 
2 r 

From ref. H we record the relation 



F' 

re = e 



ui + -r ((t> ) 



where e ^ = 1 — HlI^M. Hence, 



2{uj 



r\2 



i\2 



Arcr 



2 



27rr 



For r ~ r/i, this expression simphfies to 
-Hrh) 



2 47r ~ , / 
- + —Hrh) 



2 ■ 2 r2 



2\2 2 



where T = e^^(r) = + (0')^] 



1 



Co" I, 



1 

V 

1 



1 

4 



I? 



E_ 
V 



where 



Then the expression ( pTj ) becomes 

1 



(T;)'(r.) 



4:nV 



C 



From the field equations |^ (see Appendix A), 

V = l- 2ml 



(29) 



(30) 



(31) 



Consider for simphcity the case Vh = 1- Then, from the field equations 0] (see 
also Appendix A) 



(32) 
(33) 

(34) 
(35) 
(36) 



which is always positive because the black holes are non-extremal. (See Appendix 
A for further discussion of this point.) Thus the sign of {T^)'{rh) is the same as that 
of C. Simplifying, we have 



C = Ci + C2 + C3 + C4 + C5 



(37) 



where 



Cl = (1 - uiYuiiS - 2u;i) (38) 
C2 = ^0^(l+a;,)(l-^D(12^' + 12^^-7c^/.-9) (39) 
C3 = ^<l)l{l + uJk)\^ujl + 8uJh + 5) (40) 



C4 = -x{<Pi-vr 

1 



(41) 



C5 = ^A02,(^^-0D(l + ^/.r(2-^^ + 0D- (42) 

The first term is always positive since \uh\ < 1- The cubic in C2 possesses a local 
maximum at cuh — —0.886 where it has the value —1.724 and a local minimum at 
LVh = 0.219 where it equals —9.831. The cubic has a single root at Uh = 0.8215, and 
thus is positive for ujh > 0.8215 and negative for uj^ < 0.8215. The quadratic in C3 
is always positive and possesses no real roots. It has a minimum value of 1 when 
cuh — -1- The term C4 is always negative and C5 is always positive since \(f)h\ ^v. 

In order to assess whether or not C as a whole is positive or negative, we shall 
consider each branch of black hole solutions in turn. 

Firstly, consider the k = 1 branch of solutions. As v increases from up to 
Vmax — 0.352, Uh increases monotonically from 0.632 to 0.869. The derivative of the 
first term is 

^ = 2a;,(l - u;l)(3 - 13a;, + Sool) (43) 

where the quartic has roots given by o;^ — 1.347 or 0.278. This derivative is negative 
for LUh G (0.632, 0.869) and hence the first term decreases as v increases, and is 
bounded below by its value at the bifurcation point v — Vmax, namely 

Cl > 0.0674. (44) 

The cubic in C2 increases as Uh increases from 0.632 to 0.869, and is bounded below 
by its value when cuh = 0.632 

Uul + 12ul -7uJh-^> -5.602. (45) 

Along this branch of solutions, it is also true that 

{I-ujD < 1 - (0.632)2 = 0.601 

l+cjh < 2 

(ph < 0.19i; < 0.0669. (46) 

Altogether this gives 

C2 > -5.602 X - X (0.0669)2 ^ 3 x 0.601 = -3.767 x 10"^ (47) 
8 



The quadratic in C3 increases as uoh increases along the branch and so is bounded 
above by its value when uih = 0.869; 

AujI + 8a;fe + 5 < 14.973. (48) 

Thus, 

C3 > --^ X (0.0669)^ X 2^ X 14.973 = -7.498 x 10"^ (49) 
16 

For the fourth term, since 

{cti - v'f < (50) 

we have 

C4 > -0.15 X (0.352)^ X (^(0.601)^ + ^ x 0.15 x (0.0669)^) = -8.326 x 10~^ (51) 
and finally, 

C5 > 0. (52) 
Thus, adding these expressions up, one obtains 

C > 0.0627 > (53) 

so that (T^)'(r/i) > along the whole of the k = 1 branch of black hole solutions. 

For the quasi-A; = branch, as v decreases from Vmax down to 0, ujh increases 
monotonically from 0.869, subject to the inequality 

0.869 <ujh<l- 0.1^2. (54) 

Hence, along this branch, 

(1 - ujlf = (1 - ujhfil + Uhf > {O.lv'^y X (1.869)2 = 0.0349^;! (55) 

Thus the first term is bounded below as follows: 

ci > 0.0349^;^ x (0.869)^ x 1 = 0.0264^;^ (56) 

Since uJh > 0.869, the cubic in C2 is positive all along this branch, so that C2 and C5 
are positive. The quadratic in C3 is bounded above by its value when cuh = 2, i.e. 

Aul + 8u;h + 5< 37. (57) 

All along the quasi-A; = branch, 

(f>h < 0.5^2. (58) 

This gives 

C3 > X (0.5^2)^ X 2^ X 37 = -0.578v^ (59) 



Since v < 0.352, 



i X 0.15 X (0.5)V + (1 - (0.869)2)2 



C3 > -0.578 X (0.352)^ x = -8.874 x 10" V. (60) 
Finally, for C4 we have 

C4 > -0.l5v^ 

= -8.992 X lO^V - 2.813 x 10~V 

> -8.992 X lO^V - 2.813 x 10"^ x (0.352) V 

= -9.035 X 10" V. (61) 

In total this gives 

C > 0.0264t;^ - 8.874 x lO^V - 9.035 x lO^V = 8.491 x lO^V > 0. (62) 

In conclusion, (T^!')'(r/i) is positive for all the black hole solutions having one node 
in uj, regardless of the value of the Higgs mass v. 

Let us now check on possible contradictions with Einstein's equations. 

Consider first the case A > 1/8. Then, as r — > cxd, < and (T^)' > 0. As 
r — > r/i, < and (T^)' > 0. Hence there is no contradiction with Einstein's 
equations in this case. 

Consider now the case A < 1/8. In this case, as r 00, > and (T^)' < 0, 
whilst as r — > r^, < and (T^)' > 0. Hence, there is an interval [ra,rb] in which 
(T^)' is positive and there exists a 'critical' distance Vc G (r^, r^) at which changes 
sign. 



However, unlike the case when the gauge fields are absent |T^, here there is no 
contradiction with the result following from Einstein equations, because (T^)' > 
in some open interval close to the horizon, as we have seen above. 

In conclusion the method of ref. |T^ cannot be used to prove a 'no-scalar-hair' 
theorem for the EYMH system, as expected from the existence of the explicit solution 
of ref. . The key difference is the presence of the positive term ^ J7 in the expression 
(p!6D for {T^y. This term is dependent on the Yang-Mills field and vanishes if this field 
is absent. Thus, there is a sort of 'balancing' between the gravitational attraction 
and the non-Abelian gauge field repulsion, which is responsible for the existence of 
the classical non-trivial black- hole solution of ref. @]. However, as we shall discuss 
below, this solution is not stable against (linear) perturbations of the various field 
configurations P]. Thus, although the 'letter' of the 'no-scalar-hair' theorem of 



ref. [0, based on non-negative scalar-field-energy density, is violated, its 'spirit' is 
maintained in the sense that there exist instabilities which that the solution cannot 
be formed as a result of collapse of stable matter. 



3 Instability analysis of sphaleron sector of the 
EYMH black hole 



The black hole solutions of ref. in the EYMH system resemble the sphaleron 
solutions in SU{2) gauge theory and one would expect them to be unstable for 
topological reasons. Below we shall confirm this expectation by proving [^] the 
existence of unstable modes in the sphaleron sector of the EYMH black hole system 
(for notation and definitions see Appendix A). 



Recently, an instability proof of sphaleron solutions for arbitrary gauge groups in 
the EYM system has been given [|l^, |14| . The method consists of studying linearised 
radial perturbations around an equilibrium solution, whose detailed knowledge is not 
necessary to establish stability. The stability is examined by mapping the system 
of algebraic equations for the perturbations into a coupled system of differential 
equations of Schrodinger type |jl3|, [l^. As in the particle case of ref. the insta- 
bility of the solution is established once a bound state in the respective Schrodinger 
equations is found. The latter shows up as an imaginary frequency mode in the 
spectrum, leading to an exponentially growing mode. There is an elegant physical 
interpretation behind this analysis, which is similar to the Cooper pair instability 
of super-conductivity. The gravitational attraction balances the non-Abehan gauge 
field repulsion in the classical solution but the existence of bound states implies 
imaginary parts in the quantum ground state which lead to instabilities of the so- 
lution, in much the same way as the classical ground state in super-conductivity is 
not the absolute minimum of the free energy. 



However, this method cannot be applied directly to the black hole case, due to 
divergences occuring in some of the expressions involved. This is a result of the 
singular behaviour of the metric function at the physical space-time boundaries 
(horizon) of the black hole. 



3.1 Linearized perturbations and instabilities 



It is the purpose of this section to generalise the method of ref. [^] to incorporate 
the black hole solution of the EYMH system of ref. . By constructing appropriate 
trial linear radial perturbations, following ref. [|14|, |lT], we show the existence of 



bound states in the spectrum of the coupled Schrodinger equations, and thus the 
instability of the black hole. Detailed knowledge of the black hole solutions is not 
actually required, apart from the fact that the existence of an horizon leads to 
modifications of the trial perturbations as compared to those of ref. 
order to avoid divergences in the respective expressions 



13, 14 



m 



We start by sketching the basic steps [|T^, ^ that will lead to a study of the 
stability of a classical solution 0s (x, t) with finite energy in a (generic) classical field 



theory. One considers small perturbations 6(f){x,t) around (f)s{x,t), and specifies 
the time-dependence as 



= exp(-znt)^(x). (63) 

The linearised system (with respect to such perturbations), obtained from the equa- 
tions of motion, can be cast into a Schrodinger eigenvalue problem 

= n'^A^ (64) 

where the operators Ti, A are assumed independent of the 'frequency' Q. As we 
shall show later on, this is indeed the case of our black hole solution of the EYMH 
system. In that case it will also be shown that Ti is a self-adjoint operator with 
respect to a properly defined inner (scalar) product in the space of functions {\E'} 
[0, and the A matrix is positive definite, < \l/|y4|\l/ >> 0. A criterion for instability 
is the existence of an imaginary frequency mode in ([6^ ) 

fi^ < 0. (65) 

This is usually difficult to solve analytically in realistic models, and usually numerical 
calculations are required |jl5|. A less informative method which admits analytic 



treatment has been proposed recently in ref. |I3| , |Tl|, and we shall follow this for 
the purposes of the present work. The method consists of a variational approach 
which makes use of the following functional defined through (p3) : 



with \1/ a trial function. The lowest eigenvalue is known to provide a lower bound 
for this functional. Thus, the criterion of instability, which is equivalent to (p3|), in 
this approach reads 

fi2(^) < 
< ^1^1^ > < oo. (67) 

The first of the above conditions implies that the operator Ti. is not positive definite, 
and therefore negative eigenvalues do exist. The second condition, on the finiteness 
of the expectation value of the operator A, is required to ensure that \1/ lies in 
the Hilbert space containing the domain of Ti. In certain cases, especially in the 
black hole case, there are divergences due to singular behaviour of modes at, say, 
the horizons, which could spoil these conditions (|67|). The advantage of the above 
variational method lies in the fact that it is an easier task to choose appropriate 
trial functions that satisfy (^) than solving the original eigenvalue problem (|6^ . 
In what follows we shall apply this second method to the black hole solution of ref. 



For completeness, we first review basic formulas of the spherically symmetric black 
hole solutions of the EYMH system 0. The space-time metric takes the form [Q 



ds'^ = -N{t, r)S^{t, r)dt^ + N'^dr"^ + r'^{de'^ + sin^ Od^"^) (68) 

and we assume the following ansatz for the non-Abelian gauge potential 0, |13[ 

A = aoTrdt + aiTrdr + (cu + 1) [—r^dO + Tg sin 9dip] + ijj[rgd9 + sin 9d(p] (69) 

where u, uj and ai,i = 0, 1 are functions of t, r. The are appropriately normalised 
spherical generators of the SU(2) group in the notation of ref. [ p^ . 



The Higgs doublet assumes the form 

V2\ - #3 



^-^ft'^'^^O ; ^ = ^r (70) 



with the Higgs potential 

where v denotes the v.e.v. of $ in the non-trivial vacuum. 
The quantities u, satisfy the static field equations 



iV0" + ^0' + — 0' = -^0(^ + 1)^ + A0(0^ - t;^) (72) 
b r 2r^ 



where the prime denotes differentiation with respect to r. For later use, we also 
mention that a dot will denote differentiation with respect to t. 

If we choose a gauge in which da^ = 0, the linearised perturbation equations de- 
couple into two sectors |T3] . The first consists of the gravitational modes 6N, SS, 



6uj and 6(f) and the second of the matter perturbations 6ai, 6lj and ^ip. For our anal- 
ysis in this section it will be sufficient to concentrate on the matter perturbations, 
setting the gravitational perturbations 6N and 6S to zero, because an instability 
will show up in this sector of the theory. An instability study in the gravitational 
sector will be discussed in the following section 4. The equations for the linearised 



matter perturbations take the form |13 



(73) 

with, 

(74) 




and, 

/ Nr^ 

A=\ 2 I (75) 
V r2 

and the components of H are 

2 / , ,2 I ' i2 1 



Kiai = 2(iV5)Ma;^ + - 



= 2pl + 2NS^ — + 



4- - 

4 



Ti-a^cd = -2iNS[{p^uj) - copj^ 

n^a, = -2i[p,NSu + NS{p,uj)] (76) 

• 2 
XT 

2 2 

T-C^ai = ip*—NS(j) + i—NS{p^(j)) 
^cj-4> = ^ti>cj = —(f)NS'^ 
where the operator p^, is 

= -iNS 

dr 



p, = -iNS-^. (77) 



Upon specifying the time-dependence (|63| 



^(r,t) = ^(r)e^"* ; ^(r) = | SCd{r) \ (7J 

6ip{r) 



one arrives easily to an eigenvalue problem of the form (p^D, which can then be 
extended to the variational approach 



To this end, we choose as trial perturbations the following expressions (c.f. 

6ai = —ui'Z 

where Z is a function of r to be determined. 



One may define the inner product 



' Jr, NS 



(79) 



where is the position of the horizon of the black hole. The operator Ti is then 
symmetric with respect to this scalar product. Following ref. consider the 

expectation value 



' ' Jrn NS 



(80) 



which is clearly positive definite for real Z. Its finiteness will be examined later, and 
depends on the choice of the function Z. 

Next, we proceed to the evaluation of the expectation value of the Hamiltonian 
H (|77D; after a tedious calculation one obtains 

<^|7^|^> = / drSZ^{-2N{uj'f + 2P^N{uj^ -If 



+ PW0^(^ + 1)^ - -{oo' - ly - -<P\uj + 1)^} (81) 
4 z 

+ boundary terms 

where P = The boundary terms will be shown to vanish so we omit them in 

the expression ([8^). The final result is 



roo r 2 1 

<^|7i:|^> = / drS\-2N{u'Y--{uj^-iy--<p\u + l] 
Jr-h L r 2 



h 

+ rdr(4(^'-i)'+ 

roo f dZ^ ^ " 

Jrh V dr 



(j)\uj + iy + 2N{uj'y^S{l- Z^) 



2(^2 _ ^ ^^2^2^^ ^ 1)2 



^2) 



The first of these terms is manifestly negative. To examine the remaining two, we 
introduce the 'tortoise' co-ordinate r* defined by 

dr* 1 



dr NS 



and define a sequence of functions Zk{r*) by [ITT 



where 



Z(r*) 



k 



k = l,2, 



Z{-r*), 

1 for r* G [0, a] 



-D< ^ < 0, for r* G [a, a + 1] 



Z{r*) 



for r* > a + 1 



(83) 



(84) 



^5) 



where a, D are arbitrary positive constants. Then, for each value of k the vacuum 
expectation values of H and A are finite, < \I'|?i|\E' >< oo, and < \E'|A|\E' >< oo, 
with Z = Zk, and all boundary terms vanish. This justifies a posteriori their being 
dropped in eq. (|^). The integrands in the second and third terms of eq. (|^) are 
uniformly convergent and tend to zero as — oo. Hence, choosing k sufficiently 
large the dominant contribution in (|82D comes from the first term which is negative. 

This confirms the existence of bound states in the Schrodinger equation ([73D, (0), 
and thereby the instability ( pTD of the associated black hole solution of ref. |^ in 
the coupled EYMH system. The above analysis reveals the existence of at least one 
negative odd-parity eigenmode in the spectrum of the EYMH black hole. 



3.2 Counting sphaleron-like unstable modes in the EYMH 
system 

The exact number of such negative modes is an interesting question and we next 
proceed to investigate it. Recently, a method for determining the number of the 
sphaleron-like unstable modes has been applied by Volkov et al. [0 to the gravi- 
tating sphaleron case. We have been able to extend it to the present EYMH black 
hole. The method consists of mapping the system of linearized perturbations to a 
system of coupled Scrodinger-like equations. Counting of unstable modes is then 
equivalent to counting bound states of the quantum-mechanical analogue system. It 
is important to notice that due to the fact that the EYMH black hole solution is not 
known analytically, but only numerically, it will be necessary to make certain phys- 
ically plausible assumptions concerning certain analyticity requirements fl^ for the 



solutions of the analogue system. This is equivalent to requiring that the conditions 
for the validity of perturbation theory in ordinary quantum mechanics be applied 
to this problem. Details are described below. 

Working in the gauge Sao = 0, and denoting the derivative with respect to the 
tortoise coordinate ( ^31) by a prime, we can write the linearised perturbations in the 



sphaleron sector of the EYMH system as [13 



2N'^S^ (^uj'^ + j(f)'^^ 6ai + 2NS{u6u - u'6u) 

+ K^NS{(j)'6'^-(j)Silj') = Nr^a^6ai 
2{NSuj6aiy + 2NSuj'6ai + 5Cj" + NS'^(t)5tl) 

-1) + '^5Cj = -2a''5Cj 




]^{NSr^(P5ai)' + K^NS(l)'5ai - NS^^6u - {r^Sij')' 

+2NS' (^ ^^ + \r^\{<P' - v^)^ 6ij = r^a^dij (86) 

together with the Gauss constraint 

I (^^^i) + 2^^^ ~ ^'^'^^} = 0. (87) 

2 

Define 6^ = r6i/j, 6a = f^^ai- Then 

a^Sa = f{r*) (88) 

where 

N N r 

f{r*) = NSuj^6ai + —Sr^4>^6ai - J5Cj + u5u + —S<p5^ + - <p50 (89) 

8 4 4 



and 



f\r*) = a' [-UJ6C0 + ^05^) (90) 



and the Gauss constraint becomes 



a i6a' + u;6u - r^6n = 0. (91) 



Next, we define a 'strong' Gauss constraint by 



4 

even when a = 0. 

Using (1^) and (^2]) we may write 



5a' = -ujSCj + -<t)5^ (92) 



where 



. ±^ (%]'-^Q^Sc (93) 

Pr(f) \uj^ J Prcj) 



NS 

Pir*) = -2cu' + 2cj0' + 2cu 

r 

Q(r*) = 2^uj^ + \NS'(j)^-a\ (94) 

4 



If we substitute these expressions into the equation for 6uj we obtain the following 
equation for 5a: 



- 6a^"^ + + HP^ 5a'" 

P" 2P'2 1 
— - + 2H'P + Q-a^ + NS^J 5a" + 

\^H"P + 2P{^ +a^HP+ 



P ( I) " - (^] - i^^' + a'Q - NS'-QJ ^ fo = (95) 



where 



rP <pP 

Alternatively, we can eliminate 5E, to obtain the following pair of coupled Schrodinger 
equations p!6|: 



2. . // /20' 2iV5\ , , (2NS^ 2 AT „„\ 
(j^5a = -5a" +1-^^ j 5a + \ — —uj^ + —S^(j)^j5a + P5uj 



2 



f ^ Uuj-iy + -^\ 5u. (97) 



To proceed, it appears necessary to make the following assumptions: 

• We assume that the equilibrium solutions are continuous functions of the Higgs 
mass V. 

• We also assume that given a Schrodinger-like equation — \['" + V"^ = E"^, 
where the potential V depends continuously on some parameter v, then the 
bound state energies also depend continuously on v. This can be proven rig- 
orously if we make the physically plausible assumption of analyticity of the 
operators involved in the above system. The proof then relies on the powerful 
Kato-Rellich theorem of analytic operators fl^. This analyticity requirement 
is the case if perturbation theory is used to solve the Schrodinger system with 
potential V{v+5v) in terms of the spectrum of V{v), implying that the changes 
in the bound state energies 5E due to the infinitesimal shift in the parameter 



V are also infinitesimal. It is also the case where variational methods are appli- 
cable. This assumption implies that the eigenvalues of the discrete spectrum 
of the above equations (0), i.e. the bound-state energies for cr^ < 0, are also 
continuous as v varies continuously. 



The above assumptions had to be made because in the case of the EYMH black 
hole system, the solution is not known analytically but only numerically, and there- 
fore the issue of analyticity of the various operators involved with respect to the 
Higgs field v.e.v., v, cannot be rigorously established. 

We now notice that the continuous spectrum of the equations (0) is given by 
cr^ > 0. Hence, the number of negative modes will change by one whenever a mode 
is either absorbed into the continuum or emerges from it. 

For cr^ = the equations possess pure "gauge mode" solutions of the form 

'5« = ^, sco = -ujn, s^ = Y^ (98) 



where 



uj^ + fi. (99) 



Thus, near the event horizon, f2 ~ 0[(r — r/j)^], where /c = or 1, and at infinity, 
f2 ~ e^2 upon choosing S[oo) = 1. Hence, there is a single non-degenerate, non- 
normalisable eigenmode with = 0. 

For the fourth-order equation with = 0, 5q; ~ 0[{r — r/j)'^] near the horizon, 
where k = (twice, corresponding to the pure "gauge mode" solutions ) oi k = 

-1±V5 
2 

In the latter case, 6u,6^ ~ 0[{r — r^)^'^], and so will not remain bounded near 
the horizon. Hence there is a single non-degenerate zero mode for non-zero v. 



From the above it follows that the number of negative modes of the system (p7| ) 
cannot change at any non-zero value of v, including the bifurcation point Vmax, by 
continuity. The negative eigenvalues of this system are non-degenerate and hence 
cannot themselves bifurcate at some value of v. 



The only possible place where the number of negative modes may change is at 
V = 0. Let (p = v(f), 6C, = v6C,. Then the system of equations (P7|) becomes 



/ , 0' NS\ , 
\ 6 r 



-5Cj" - 2 



4 



5uj 
5a 



a 5a 



a^5Cj (100) 



together with the Gauss constraint 5a' = —uj5uj + jv'^r(j)5^. 



Consider the k = n branch of the EYMH solutions. In this branch, ^ has a 
well-defined limit as f ^ 0, and the system (^) is continuous at w = 0. In this case, 
the Gauss constraint reduces to 5a' = —uj5ijJ , and substituting in ( |10U| ) yields the 
equation 



9 9 A^^2 

5a H — UJ 5a H — uj 5a = a 5a. 

UJ 



'101] 



This is the equation studied in ref. where it was shown that there are exactly n 
negative eigenvalues. Furthermore, there is a single non-degenerate zero mode given 
by 

r^Q' , , 

5a = TTTTT^ (102) 



2NS^ 



where 



/rHV_ 
[2NS^ 



uj^VL. 



(103) 



As before, near the horizon ~ 0[(r — r/j)] or 1] ~ 0[1] whilst at infinity i7 ~ r 
or r2, giving a single eigenmode with zero eigenvalue. Thus, the number of negative 
eigenvalues does not change at f = for this branch of solutions. 



For the quasi k = n — 1 branch of solutions 4- does not have a well-defined limit 
as u — > 0, and the system (^) is not continuous at f = 0. Hence, by continuity, we 
can conclude that both the k = n and the quasi k = n — 1 black holes have exactly 
n unstable modes in the sphaleron sector. 



4 Instabilities in the Gravitational Sector - 
Catastrophe theory approach 

It is the aim of this section to prove the existence and count the exact number of 
unstable modes in the gravitational sector of the solutions. In the first part we shall 
study the conditions for the existence of unstable modes in a linearized framework, 
and we shall study the possibility of a change in the stability of the system as 
one varies the Higgs v.e.v. v continuously from up to the bifurcation point (c.f. 



figure 2). In the second part, which will deal with the change of the stability of 
the system at the bifurcation point, we shall go beyond linearized perturbations by 
applying catastrophe theory. It should be stressed that although catastrophe theory 
was first employed by the authors of 0, our approach in this section is somewhat 
different, and has certain advantages, not least of which is that we are able to exploit 
the known stability of the Schwarzschild black hole to draw conclusions about the 
non-trivial EYMH black holes. 



4.1 Linearized perturbations 

The linearised perturbation equations for the gravitational sector are: 



2NS 

b 0(p +Ud,uOuJ + 

r 



(104) 



where the prime denotes differentiation with respect to the tortoise coordinate (|83|) 
as before, and 



U.., 









3uj 













1 + -r - 4r2cj'2 
4 



{NSy\ uj{uj^-l) , 



r ^ S 



+ 



-ijj + 2ruj 



'N^iNSn 

T SI 



+ 2r(f)'uj{uj^ - 1 



Us 



, r S 



+20X1 + uyr + AX(f)(j)'r\(f)^ - d^) 



(105) 



The continuity argument described previously when applied here, implies that the 
number of negative eigenvalues of the system ( |104| , p.05| ) can change only when there 
is a zero mode. 



Suppose that for some v there is such a zero mode of ( p.04|) . Given the 
background solution (0^,0) of the field equations, then (tu + 6uj,(f) + 6(j)), together 
with the corresponding metric functions, will also be a solution of the field equations. 



As r oo , then 
and 



- 0[e-'"'*] r^-oo (106) 
~ const r —>■ — oo. (107) 



Also, since in this sector the change in the mass function m(r) is given by (see 
Appendix A) 

6m(r) = 2N^6uj(r) + r^N^6(l)(r) (108) 
dr dr 

we have that for this zero mode 6m — as r* — ±oo. Hence, for fixed param- 
eters rh,X,v,g there are two solutions of the original field equations ^ ^ ^5 ^'^d 

^ ^ ^0 ^ ) satisfying the required boundary conditions. For our purposes we shall 

assume that the solution of ref. is unique. Compatibility with the above analy- 
sis, then, requires the latter to be valid only at the bifurcation point (c.f. figure 2). 
Hence, there is a single zero mode at f = Vmax- For any other v the zero mode is 
absent. 



Let 6(f> = vS(j), (j) = v(j). The equations ( |1(J5| ) become 

f + U^Jco + U^^ScP = aH(t) (109) 



r 

where U^^ = vU^^, If^^ = vU^^. 

These equations have a well-defined limit as f — > 0, and are continuous at f = 0. 
At i; = the equations reduce to 

Soo" + U^Juj = a'^Suj (110) 

where 



.,,N , iNSy\ , uiuj'-l)^^, 



3uj'-l- Arico'Y ( - + 1 + ^ 'Su' 



If this equation possesses a zero mode, then 600 —>■ const, as r — >■ 00 for this mode. 
Asr ^ 00, ^ l-^ + 0[e-T, 00 -l + 0[e-% S l + Ofe^, and the equation 
takes the form 

^(*-) = -;3 1-- 5;**-) + ;i ( 1 - - ) ("2) 



Suppose that 6uj = J2'^=o(^nr oq 7^ as r — > 00. Then p = 2 or —1. Let 
6uj = r'^f{r). Then, /(r) satisfies the equation 

r^f(i-'^y(i-'-^)^o. (113) 



The solution to this equation as r ^ 00 assumes the form 

fr — 1M.\ A A 



with A, B arbitrary constants and Ai 



linir^oo '^('")! so that 



Ar'^ 



log 



( 



r-2M 



) 



Ar 



A 



(115) 



+ 



+ 



AM 



r 



Hence 6uj can remain bounded as r ^ oo only if A = i? = 0, i.e. only the trivial 
solution exists. Thus, there are no non-trivial zero modes for w = 0. 

We, therefore, conclude that along each branch of solutions the number of negative 
modes remains constant from f = to f = Vmax- 

4.2 Bifurcation points and catastrophe theory 

In order to determine what happens at f = Vmax "we appeal to catastrophe the- 
ory . Our aim is to study the possibihty of a change of the stability of the system 
at Vmax- To this end, we have to determine a certain function {catastrophe func- 
tional) in the black hole solution which changes dis continuously despite the smooth 
change of certain (control) parameters of the system. As we shall show below, in the 
case at hand the role of the catastrophe functional is played by the mass function of 
the black hole, whilst the control parameter is the Higgs v.e.v. v. At the bifurcation 
point Vmax we shall find a fold catastrophe which affects the relevant stability of 
the branches of the solution. In addition, the catastrophe-theoretic approach allows 
for an exact counting of the unstable modes in the various branches. For notations 
and mathematical definitions on catastrophe theory we refer the interested reader 
to Appendix B. 

We should note at this point that although catastrophe theory seems power- 
ful enough to yield a universal stability study of all kinds of non-Abelian black 
holes 0, however one should express some caution in drawing conclusions about 
absolute stability. Indeed, catastrophe theory gives information about instabilities 
of certain modes of the system. If catastrophe theory gives a stable branch of so- 
lution, this does not mean that the system is completely stable, given that there 
may be other instabilities in sectors where catastrophe theory does not apply. In 
our EYMH system this is precisely the case with the sphaleron sector. However safe 
conclusions can be reached, within the framework of catastrophe theory, regarding 
relative stability of branches of solutions, and it is in this sense that we shall use 
it here in order to count the number of unstable modes of the various branches of 
the EYMH system. Having expressed these cautionary remarks we are now ready 
to proceed with our catastrophe-theoretic analysis. 

The mass functional M. (c.f. appendix A) can be re-written as a functional of the 
matter fields only as follows: 



First note that Ai = m(oo). Let /i(r) = m(r) — m{rh) = m{r) — -y. Then, using 
a prime to denote d/dr 



^'{r) = m{r) 



2m 



'\2 I ^2/i/\2\ 



r 



1 (2(w') + r 

r / 

^ ^ + -^(1 + + (02 _ ^2)2 



rh 



2r' 



I „2/ j/\2\ 



+ - 



+ _(02_~2)2 



r 



:ii6) 



The last term on the right-hand-side can be written in terms of the metric function 
6 (cf appendix A) 



-l^{2{u'y + r 
r 



t\2 I „2lit\2\ 



Solving for /x gives 



(117) 
(118) 



where 



}C[uj,(t>] = 



{2{ur+A<p'r 



\ (1 — uj 



2\2 



+ 



2r2 



Hence, setting 6{oo) = we obtain 

M = —+ r }C[iu,^]e-^^''Ur. 

2 Jru 



(119) 



(120) 



Varying this functional with respect to the matter fields yields the correct equa- 
tions of motion [|18]. Thus, the equilibrium solutions of the field equations will be 
stationary points of the functional M.. 



If we plot the solution curve in {y, 5q, M.) space, then the resulting curve is smooth 
(cf. figure 1). For the Catastrophe Theory (cf. Appendix B) we consider 5q as a 
variable, and f as a control parameter. The Whitney surface could be defined in 
our case as follows. For each f , consider a smoothly varying set of functions cj^, 0^ 
indexed by the value of 5q they give: 

/•CO 1 

5o= / -{2{u',f + r\<P',f)dr (121) 

such that a;^, 4>s are the appropriate solutions to the field equations when 6o lies on 
the solution curve. Then, the solution curve represents the curve of extremal points 
of this Whitney surface, SAi = 0. 



The projections of this curve onto the {v,6q) and {60, Ai) planes are also smooth 
curves. The catastrophe map x projects the solution curve onto the (f, plane 

X : {v,So,M)^iv,M). (122) 

This yields the curve shown in figure 2. This map is regular except at the point 
(f = Vmax,-M. = A4max), whcrc it is siugular. This point is the bifurcation set B. 



Since the Whitney surface describes a one-parameter [v) family of functions of 
a single variable (6), and the bifurcation point is a single point, we have a fold 
catastrophe, as found in ref. from a different point of view. A more detailed 
comparison, of the results of that reference with ours will be made at the end of the 
section. 



Catastrophe theory tells us that the stability of the system will change at the point 
B on the solution curve, and, furthermore, that the branch of solutions (including 
the point B) having the higher mass (for the same value of v) will be unstable, 
relative to the other branch. Hence, from our previous continuity considerations, 
the k = n branch of solutions will have exactly one more negative mode than the 
quasi- k = n — 1 branch. 



The catastrophe theory analysis applies to the gravitational sector rather than the 
sphaleronic sector, since gravitational perturbations correspond to small changes in 
the functions u and 0, whilst keeping the functional form of A4 fixed. On the 
contrary, sphaleronic perturbations keep the functions u and fixed, affecting the 
functional form of A^. As we discussed in the previous section, the number of 
unstable sphaleron modes is the same for the k = n and quasi-fc = n — 1 branches. 



All that remains is to determine the number of negative modes of the quasi-A; = 
branch of solutions. From the above considerations, this will be equal to the number 
of negative modes of the v = limiting case of this branch of solutions, which is 
nothing other than the Schwarzschild black hole. The gravitational perturbation 
equation is in this case, where a prime now denotes d/dr*, 



5uj" + U^ujSuj = a'^duj 



(123) 



where 



\r ^ J f 



(124) 



For the Schwarzschild solution A^ = 1 — ^,5* = 1,uj = 1. Equation ( |123| ), then, 
reduces to 

- 5tu" + ^ (^1 - ^) 5cu = a''5uj (125) 



which has the form of a standard one-dimensional Schrodinger equation with poten- 
tial 

V r / ar* \ r J dr 
As r* ^ ±cxD, ^(r*) — >■ 0. On the other hand, for finite r*, —oo < r* < oo, the 
potential is positive definite, V{r*) > 0. 



Then, by a standard theorem of quantum mechanics , the Schrodinger equation 
( |125 , 126| ) has no bound states. Thus, the Schwarzschild solution (and, hence, the 
quasi-/c = branch of solutions) has no negative gravitational modes, a known result 
in agreement with the no-hair conjecture. 



Working inductively through the various branches of solutions (the v = limit of 
the k = n branch is the same as that for the quasi-/c = n — 1 branch after replacing 
uj by —uj) we find that the k = n branch possesses exactly n unstable gravitational 
modes, and the quasi-A; = n — 1 branch exactly n — 1 negative modes. This result 
has been conjectured but not proven in ref. 

Before closing the section we would like to compare our results with those of ref. 
p. In ref. the authors also used catastrophe theory to draw conclusions about 
the stability of the EYMH black holes, but their approach was somewhat different. 
There the writers fixed the parameter A = 0.125 and also fixed the Higgs mass v. 
They then varied the horizon radius Vh and for each solution calculated the value of 
the mass functional Ai, the field strength at the horizon Bh given by 



horizon 

(127) 



and the Hawking-Bekenstein entropy S = TiVh, to give a smooth solution curve in 
{M.,Bh,S) space. The projection of this curve on to the {M.,S) plane has the 
same qualitative features as our figure 2. Here we have also fixed A = 0.15, and 
in addition we have fixed r/^ = 1 and varied the Higgs mass v from up to the 
bifurcation point. Torii et al concluded that the k = 1 branch of solutions was 
more unstable than the quasi-fc = branch of solutions. The advantage of our 
approach is that, by interpolating between the various coloured black hole solutions 
1l9[] , beginning with the Schwarzschild solution, we will be able to calculate the exact 
number of unstable modes of each branch of solutions and not just give qualitative 
information concerning their relative stability. 



5 Entropy considerations 

It remains now to associate the above catastrophic considerations with some ele- 
mentary 'thermodynamic' properties of the black hole solutions, and in particular 
their entropy. Our aim in this section is to give elementary estimates of the entropy 
of the various branches, assuming thermodynamic equilibrium of the black hole with 



a surrounding heat bath. Such estimates will allow an association of the stability 
issues with the amount of entropy carried by the various branches of the solution. In 
particular we shall argue that the 'high-entropy' branch has relatively fewer unsta- 
ble modes than the 'low-energy' ones and, thus, is relatively more stable. We shall 
employ approximate WKB semi-classical methods for the evaluation of the entropy. 
We shall also study the conditions under which such estimates are valid. 

The calculation of the entropy of the black hole will be made on the basis of 
calculating the entropy of quantum matter fields in the black hole space time. This 
will constitute only a partial contribution to the total black hole entropy. A com- 
plete calculation requires a proper quantization of the gravitational field, which at 
present is not possible, given the non-renormalizability of (local) quantum gravity. 
Ignoring such back reaction effects of the matter fields to the (quantum) geome- 
try of space time results in ultra-violet divergences in the calculated entropy of the 
matter fields J^, If. Such divergences can be absorbed in a renormalization of the 
gravitational (Newton's) constant. This is so because the entropy is proportional 
to the area of the black hole horizon, with the divergent contributions appearing as 
multiplicative factors. 

In what follows we shall estimate the entropy of a scalar field propagating in the 
EYMH black hole background. Anticipating a path integral formalism for quantum 
gravity, we shall compute only the entropy which is due to quantum fluctuations of 
the scalar field in the black hole background. The part that contains the classical 
solutions to the equations of motion contributes to the 'classical' entropy, associated 
with the classical geometry of space time. This part is known to be proportional 
to 1/4 of the horizon area 0, ||. The quantum-scalar- field entropy part will also 
turn out to be proportional to the horizon area, but the proportionality coefficient 
is linearly divergent as the ultraviolet cut-off is removed, exactly as it happens 
in the corresponding computation for the Schwarzschild black hole [0. Absorbing 
the divergence into a conjectured renormalization of the gravitational constant will 
enable us to estimate the entropy of the various branches of the EYMH black hole 
solution, and relate this to the above-mentioned catastrophe-theoretic arguments. 

As we shall show, this will be possible only in the non-extremal case, which is the 
case of the numerical solutions studied in ref. |^ and in the present work. Among 
the solutions, however, there exist some extremal cases, for which the Hawking 
temperature - which is defined by assuming thermal equilibrium of the black hole 
system with a surrounding heat bath - vanishes. In such a case, the linearly divergent 
entropy of the scalar field vanishes. However, there are non trivial logarithmically 
divergent contributions to the black hole entropy which cannot be absorbed in a 
renormalization of the gravitational constant. Moreover, the classical Bekenstein- 
Hawking entropy formula seems to be violated by such contributions to the black- 
hole entropy. The situation is similar to the case of a scalar field in an extreme 



(3+l)-dimensional dilatonic black hole background |p, and seems to be generic to 
black holes with non-conventional hair. We shall briefly comment on this issue at 
the end of the section. 



We shall be brief in our discussion and concentrate only in basic new results, 
relevant for our discussion above. For details in the formalism we refer the interested 
reader to the existing literature 0, To start with, we note that the metric for 
the EYMH black holes is given by: 



ds^ = - 1- 



2m(r) 



e-^'(^)dt^ + 1 - 



2m(r) 



-1 



dr^ + r\d9^ + sm^9d^^). (128) 



Consider a scalar field of mass propagating in this spacetime 0, satisfying the 
Klein-Gordon equation: 



-9 



(129) 



Since the metric is spherically symmetric, consider solutions of the wave equation 
of the form 

^t,r,e,ip)=e-''''fEiir)Yi^,ie,cp) (130) 

where Yim^{6, (p) is a spherical harmonic and E is the energy of the wave. The wave 
equation separates to give the following radial equation for fEiij") 



1 



2m(r) 



3^d_ 
J.2 ^j. 



-Sir)^2 A _ 2m(r)\ dfEijr) 
\ r I dr 



/([+!) 

2 



fEiir) 



[131) 



The "brick wall" boundary condition is assumed 0, namely, the wave function is 
cut off just outside the horizon. 



$ = at r = r/t + e 



(132) 



where is the black hole horizon radius, and e is a small, positive, fixed distance 
which will play the role of an ultraviolet cut-off. We also impose an infra-red cut-off 
at a very large distance L from the horizon: 



$ = at r = L, where L ^ r^. 



(133) 



Hence / satisfies 

fEi{f) = when r = + e oi r = L. 



(134) 



In anticipation of being able to use a WKB approximation, define functions K{r) 
and h{r) by 



h{r) 



2m(r) 



-1 



E' 1 



2m(r) 



1(1 + I) 2 
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2m(r) 



(135) 
(136) 



Then the equation for fEiij) becomes 



1 d 



h{r) dr 

Now define a function u{r) by 



fEi{r) 



+ K\r,l,E)fEi{r) = 0. 



■u(r) 



h{r) 



Then u(r) satisfies 



dr"^ 



+ 



1 fdh" 
Ah^ \ dr 



1 d'^h 
2hd^ 



u = 0. 



The WKB approximation for u will be valid if 



1 fdh' 



1 d'^h 



2hdr'^ 



and 
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(137) 
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(139) 



(140) 



The first inequality is required so that u can be taken to satisfy the equation 

d^U ^^n 

— + = 



:i4ii 



(142) 



where K is now the radial wave number, and the second inequality is required so 
that the approximation to the wave function 



u{r) 



: exp 



±i J K{r)dr 



(143) 



is valid. Assuming, for the present, that the WKB approximation is valid, define 
the radial wave-number K as above whenever the right- hand- side of ( |135| ) is non- 
negative. Define = otherwise. Then the number of radial modes tik is given 
by 



nriK 



drK{r, I, E) 



(144) 



where the fact that uk must be an integer restricts the possible values of E. For 
fixed energy E, the total number N of solutions with energy less than or equal to 
E is 



where the integration is performed over all values of I such that the argument of the 
square root is positive. 

The Hawking temperature of the black hole is given by 

T-'=P=^^ (146) 

where 6h is fixed by the requirement that 6{oo) = 0, and ' = d/dr. Assume further 
that <^ 1. It should be noted that T = in the extreme case m'^j = 0.5. Further 
comments on the entropy in this situation will be made at the end of the section. 
However, as discussed in Appendix A, this situation does not arise for the black 
holes we are concerned with. The free energy F of the system is given by 





(145) 
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(147) 
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1 - ex.-p(-(3E) ■ 



Hence 



log(l - e-^^) 



nK,l,mi 




(148) 



for large (3, integrating over appropriate I, E. Integrating by parts, 



F 




(149) 



where we have substituted for hk from ( |144| ). The / integration can be performed 
exphcitly, 



dl{2l + 1) 
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to give 



F = — 



dE 



Svr J exp{(3E) — 1 Jr^ 



drr 1 — 



2m(r) 
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r 



E' 



' ^ 2m(r) 



Introduce a dimensionless radial co-ordinate x by 

r 



a; 



51) 
(152) 



Then 
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dE- 



Stt j exp{f3E) — 1 Ji+e 
where e = — , L = — , and mfx) = "''-^''^^ . 
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fx 



(153) 



The contribution to F for large values of x is 



Fo = 

971 



dE 



(154) 



exp{(3E) - 1 

which is the expression for the free energy in flat space. The contribution for x near 
1 diverges as e — > 0. For x near 1, the leading order term in the integrand in (|153|) 
is 

E^{x-l)-\l-2mf^f (155) 
where rh'f^ = tti' (1) = m' (rh) . This gives the leading order divergence in F 
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dE 
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The total energy U and entropy S are given by 

^ = dp^^^^ = 157 p-^ 



(156) 
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S = f3 



.OF 8 



Substituting for /3 from (|4 



dp Abe 
obtain 



/33 



S 



360e" 



[l-2m')e 



-3<5h 



(158) 



(159) 



Before discussing the implications of this formula, it is necessary to ascertain 
when the approximations used are valid. Firstly, /5 ^ 1 if r/^ ^ 1 or 1 — 2m'^ <C 1. 
In the first case, non-trivial (viz. non-Schwarzschild) solutions exist only for very 
small values of f , the Higgs mass and these solutions will be very close to the 
Schwarzschild solution having mass M = rh/2. In the second case, the black hole 
is very nearly extremal. For 1 = 2m'f^ exactly, the above analysis does not apply. 
However, 1 — 2m'^ -C 1 for large n [19| and for any value of v for which a non-trivial 



solution exists. We shall discuss the physical implications of this (nearly) extremal 
case at the end of the section. 

Secondly, consider the validity of the WKB approximation. The principal contri- 
bution to the free energy F comes from the region where K is large; in particular, 
above we have concentrated on x close to 1. It is expected that the WKB ap- 
proximation will be valid when K is large. For K large, it may be approximated 
by 

2m(r)\ 



For r near r^, then 



whence 



K = E{1 

\ r J 

K = E{1- 2m;)-^(r - Th)'^ + 0(1) 
dK 



Hence 
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< \K^\ if 



E 



1 - 2mi 



> 1. 



Similarly, for r near we may approximate h by 

h = e-V(l - 2m;) (r - Vh) + 0(r - Vhf. 

Then 
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so that 
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(160) 

(161) 
(162) 
(163) 

(164) 
(165) 
(166) 



Thus, for black hole solutions with large n, the WKB approximation is valid except 
for small values of E. Now return to the expression for the entropy (|158|) , 



S = Sunear = ^(1 " 2m'^)e-^'\ (167) 

We notice first that the entropy is positive, due to the fact that for the solutions 
m'^ < 1/2 to avoid naked singularities. Having said that, we now fix n and consider 
the two branches of black hole solutions, the k = n and quasi-fc = n — 1 solutions. 
The linear divergence r^j/e is a common multiplicative factor in all branches, and 
thus can be absorbed in a renormalization of the gravitational constant . This can 
be done as follows: Re-write r/i/e = 47rr^^^^, where A = 47rr^ is the horizon area. 
Go is the bare graviational coupling constant (which, by convention, had been set to 
one in the previous formulae), and r^e = 2mhe may be considered as the invariant 
distance (cut-off) of the brick wall from the horizon. The classical Bekenstein- 
Hawking entropy formula is then still valid, but with the renormalized gravitational 
constant Gr replacing the bare (classical) one Go 

Sclassical + Sunear = {~rFr + ^[~])^ = ^ (168) 

4(jro e 4(jr^ 



Such a renormalization may be thought of as expressing quantum matter back 
reaction effects to the space-time geometry. Doing this in our case, we observe 
from ( |167|) that for each f , the k = n solution has larger m')^ and 6h than the 
quasi-A; = n — 1 solution. Hence the k = n solution has a lower entropy than the 
quasi- /c = n — 1 solution, in agreement with Torii et al [0]. 



Before closing the section we would like to make some important comments con- 
cerning the extreme case m'{rh) = 1/2, for which the Hawking temperature (|146|) 
vanishes. In this case the linearly- divergent part of the entropy (|159| ) also vanishes, 
but this is not the case for the next-to-leading order logarithmically divergent partQ. 

The logarithmic divergent part of the free energy can be found from ( |153|) by 
requiring the following expansion 

m(x) = rhh + rh'hix — 1) + -m'/^(x — 1)^ + . . . rhh = -. (169) 

Using the trick x = 1 -|- (x — 1) we can write down the identities: 

x"^ = l-(x-l) + (x-l)^ + ... 

= l + 2(x-l) + (x-l)^ (170) 

^It should be noted that the logarithmic divergent parts exist also in the non-extreme case, but 
there they are suppressed by the dominant linearly divergent terms. It can be easily checked that 
for the solutions of ref. their presence does not affect the entropy considerations above, based 
on the linearly divergent term. 
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(171) 



Substituting in ( |153| ) we obtain for the next-to-leading order divergence of the free 
energy 
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(172) 



Using the formulae 
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(173) 



the expression (|172|) reduces to 



4 ,7r4 2-4mi+m'/ 
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(174) 



From 



the corresponding next-to-leading contribution to the entropy in the 



extremal case m', = 1/2 (where the linear divergence vanishes) is given by the 
following expression : 



Snlo 



Uli- 



log — 



(175) 



Thus, we observe that in the extremal case the entropy diverges logarithmically 
with the ultraviolet cut-off, in a similar spirit to the case of the dilatonic black hole 
background In our case, however, the horizon area does not vanish, because 
there is no dilaton field exponentially coupled to the graviton. Thus, one could hope 
that the divergent contribution ( |175| ) could be absorbed in the renormalization of 
the gravitational constant, so that a formal Bekenstein-Hawking expression for the 
entropy is still valid. However, as we see from (|175|) , for generic scalar fields this 
cannot be the case, due to terms that spoil the proportionality of Snio to the black 
hole horizon area A = ^irr^. Indeed, let us analyze the various contributions in 



The first term can be absorbed into a renormalization of tlie gravitational con- 
stant, and respects tlie classical formula ( |168|) . This is not the case with the second 
term however. From equation ( |202D of Appendix A, we observe that there are con- 
tributions that depend on the (boundary) horizon values of the fields (ph and ujh 
which are not proportional to the horizon area A, 

ml = -(l-^{l + u;,r] (176) 

Thus, the associated contribution to the black hole entropy seems not to be related 
to geometric aspects of the black hole background. 



One is tempted to interpret such contributions as being associated with infor- 
mation loss across the horizon. This is supported by the fact that the logarithmic 
divergencies disappear for black holes whose horizon is vanishing in the sense of 
rh e. For consistency with the interpretation as loss of information, the positivity 
requirement of the relevant contribution to the black hole entropy has to be imposed. 
Returning to formula ( |175| ) we observe that the above requirement implies m'l > 0. 
In the present case we do not know whether extremal solutions exist. A priori there 
is no reason why such solutions should not exist in the EYMH system. If such a 
solution exists, the above-mentioned positivity requirement will impose restrictions 
on the boundary (horizon) values of the hair fields of the black hole background. 
From the case at hand, it seems that the ambiguities in sign are associated with 
the presence of the non-abelian gauge field component uj. Indeed, from (|176|) , it is 
immediately seen that the contributions of the scalar Higgs field alone to m'l^ are 
manifestly positive. The terms that could lead to negative logarithmic contributions 
to the entropy are associated with the field uj and vanish for a; = —1. 



This phenomenon is somewhat similar to what is happening in the case of a 
(spin one) gauge field in the presence of an ordinary black hole background. If one 
integrates quantum fluctuations of a spin one field in a gravity background, there 
is an induced coefficient in front of the Einstein curvature term in the effective 



action whose sign is negative for space-time dimensions less than 8 ||22[. Notice 
that such sign ambiguities do not occur for scalar fields in conventional black hole 
backgrounds. In our case, there are gauge fields present in the black hole background 
associated with non-conventional hair. The sign ambiguities found above in the 



logarithmically-divergent contributions to the entropy (|1 75|) occur already when 
one considers quantum fluctuations of scalar fields. This is associated with negative 
signatures of terms that involve the gauge field hair background in the effective 
action. 



Some comments are now in order concerning the the so-called entanglement en- 
tropy [^] of fields in background space times with event horizons or other space-time 
boundaries. The entanglement entropy is obtained from the density matrix of the 



field upon tracing over degrees of freedom that cross the event horizon or he in the 
interior of the black hole, and therefore is closely associated with loss of information. 
The entanglement entropy is always positive. This immediately implies a difference 
from the ordinary black hole entropy, computed above, for the case of spin one 
fields [^. On the other hand, for scalar fields in ordinary black hole backgrounds 



both entropies are identical, since in that case sign ambiguities in the entropy do not 
arise. On the other hand, our computation for the extreme EYMH case, provided 
the latter exists, has shown that, in general, one should expect a difference between 
the two entropies even in the case of scalar fields propagating in such (extreme) 
non-Abelian black hole backgrounds. 

There exists, of course, the interesting possibility that the entanglement entropy 
of scalar fields in this extreme black hole background can be identified with the 
logarithmic entropy terms ( |175| ), in which case the latter must be positive definite. 



This, as we discussed above, would imply restrictions on the boundary (horizon) 
values of the gauge hair for the extreme black hole to exist. The restrictions seem to 
be relatively mild though. As an example of the kind of the situation one encounters 
in such cases, consider the case where extreme EYMH black hole solutions exist. 
From (|176| ), we observe that positivity of m'l^ implies restrictions on the size of ujh, 
2ujh > + > 0, which is a mild restriction. 

However, all these are mere speculations at this stage. One has to await for 
a complete analytic solution of the EYMH black hole problem before reaches any 
conclusions regarding entropy production and information loss in extreme cases. 
Therefore, we leave any further considerations on such issues for future work. 



6 Conclusions and Outlook 

In this work, we have analyzed in detail black holes in (3+l)-dimensional Einstein- 
Yang-Mills-Higgs systems. We have argued that the conditions for the no-hair theo- 
rem are violated, which allows for the existence of Higgs and non-Abelian hair. This 
analytic work supports the numerical evidence for the existence of hair found in |^ . 
This is due to a balance between the gauge field repulsion and the gravitational 
attraction. However we have shown that the above black holes are unstable, and 
therefore cannot be formed by gravitational collapse of stable matter. Although 
the instability of the black hole sphaleron sector was expected for topological rea- 
sons, however our analysis in this work, which includes an exact counting of the 
unstable modes in this sector, acquires value in the sense that we have managed to 
describe rigorously the sphaleron black holes from a mathematical point of view. In 
the gravitational sector we have used catastrophe theory to classify and count the 
unstable modes. Our method of using as a catastrophe functional the black hole 
mass and as a control parameter the Higgs field v.e.v. proved advantageous over 



existing methods of similar origin in that we managed to understand the connec- 
tion with the Schwarzschild black holes from a stability/catastrophe-theoretic point 
of view. The above analysis, although applied to a specific class of systems, how- 
ever is quite general and the various steps can be applied to other self-gravitating 
structures in order to reach conclusions related to the existence of non-trivial hair 
and their stability. For instance, we can tackle the problem of moduli hair of black 
holes in string- inspired dilaton-coupled higher derivative gravity |24]. The presence 
of Gauss-Bonnet combinations in such systems shares many similarities with the 
case of the non-Abelian black holes, and it would be interesting to study in detail 
the possibility of having non-trivial hair (to all orders in the Regge slope a') and its 
stability, following the methods advocated in the present work. 



In addition to the question of the stability of non-conventional hairy solutions, 
the above analysis has revealed another important aspect concerning the information 
theoretic content of these (3+l)-dimensional hairy black holes, namely the existence 
of logarithmic divergent contributions to the entropy of matter (quantum (scalar) 
fields) near the horizon. Such contributions owe their existence to the non-trivial 
hair of the black hole, and they modify the Bekenstein-Hawking entropy formula, 
by yielding contributions that do not depend on the horizon area. Our findings can 
be compared to a similar situation characterizing extreme (string-inspired) black 
holes 0. There, the deviation from the Bekenstein-Hawking entropy was seen to 
occur by the fact that in the extreme case, due to the presence of the dilaton, 
the effective horizon area vanishes, whilst the entropy did not vanish. In our case, 
despite the non vanishing entropy in the extreme case, the logarithmically-divergent 
entropy contributions violate explicitly the classical entropy-area formula by yielding 
contributions that are independent of the horizon area. This kind of entropy is 
clearly associated with loss of information across the horizon but it is not described 
in terms of classical geometric characteristics of the black hole. If true in a full 
quantum theory of gravity, this phenomenon might explain the information paradox. 
The question of associating this entropy with the entanglement entropy of fields in 
the EYMH background is left open in the present work. We hope to come back to 
this issue in the near future. 



Whether a full quantum theory of gravity could make sense of such divergencies 
or not remains to be seen. There are conjectures/indications that string theory, 
which is believed to be a mathematically consistent, finite theory of quantum grav- 
ity, yields finite extensive quantities at the horizon [§, if string states, which in a 
generalized sense are gauged states, are properly taken into account |2^. However, 
our understanding of these issues, which are associated with the incompatibility - 
at present at least - of canonical quantum gravity with quantum mechanics, is so 
incomplete that any claim or attempt to relate the above issues to realistic com- 
putations involving quantum black hole physics would be inappropriate. We think, 
however, that it is interesting to point out yet another contradiction of quantum me- 
chanics and general relativity associated with the proper quantization of extended 



objects possessing space-time boundaries. 
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Appendix A 

Notation and conventions 

Throughout this paper we use the sign conventions of Misner, Thorne and Wheeler 
[p^ for the metric and curvature tensors. In particular, the signature of the metric is 
( — h ++). For the EYMH system, we write the most general spherically symmetric 
metric in the form 

ds^ = -NS^df + N-^dr^ + r\de^ + sin^ Od^^) (177) 

where and S are functions of t and r only and can be written in terms of the 
mass function m and the function 6 as 

iV(t,r) = l-^^:^^, Sit,r) = e-'^'''-l (178) 
r 

This latter form of the metric is particularly useful for black hole space-times. Fol- 
lowing ref. [^, we take the most general spherically symmetric SU(2) gauge potential 
in the form 

A = aoTrdt + aiTrdr + (1 + cu) [rg sin 9d(p — r^dO] + Cj[TQd6 + sin Odip] (179) 

where Oq, cti, and Cj are functions of t and r alone and the are given by 

Tr = Ti sin 6' COS + r2 sin 6' sin + COS 6* (180) 
Tg = Ti COS COS <y9 + r2 COS sin (y9 — Ts sin ^ (181) 
= —Ti sin ip + T2 cos (f (182) 

with Ti, i = 1,2,3 the usual Pauh spin matrices. The complex Higgs doublet assumes 
the form 

$ = ^f^;+.^;^0 (183) 

where a suitable spherically symmetric ansatz is 

ip = ilj{t,r)r, (f) = (f){t,r). (184) 
Then the EYMH Lagrangian is M 



^EYMH — — — 

An 



:i85) 



and the Higgs potential is 

V{<l>') = ^{<P'-vr. (186) 



For the equilibrium static solutions, and ijj all vanish and the remaining 

functions depend on r only. The metric functions m(r) and 6{r) are required by the 
Einstein equations to satisfy the following, where ' = d/dr: 
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2 in order for a regular event horizon 



subject to the boundary conditions m{rh) ~ — 
at r = r/j, and, in order for the spacetime to be asymptotically flat, 5(oo) = 0. For 
an asymptotically fiat spacetime, it is also the case that m(r) — M as r — oo, 
where M is a constant equal to the ADM mass of the black hole. Integrating (|187|) 
from Th to oo we obtain: 

y = m(oo) - m{rh) -- 
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This equation defines the mass functional Ai as an integral of the fields over the 
spacetime. 



Finally we define the 'tortoise' co-ordinate r* by 

dr* _ 1 



(190) 



Numerical solution of equilibrium equations 



The static field equations for the metric functions and matter fields are: 
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(193) 
(194) 



For finite energy solutions, we require that uj{oo) = —1, (pioo) = v and 6{oo) = in 
order for spacetime to be asymptotically flat. These equations trivially possess the 
Reissner-Nordstrom solution given by 

m = y, uj = -l, (p = v, 6 = 0. (195) 

Non-trivial solutions do not occur in closed form, so a numerical method of solution 
is necessary as in Q]. We set the horizon radius rh = I and A = 0.15 (cf. A = 0.125 
in ref. |[). 

From the above equations, if the function 6{r) satisfies ( |192| ) then 6{r) + constant 
will also be a valid solution. To make the numerical solution easier, we set S{rh) = 
(so that S{oo) = will not be satisfied) when integrating outwards from r/j. An 
appropriate constant can then be added to 5(r), after the field equations have been 
solved, to ensure that the boundary condition at infinity holds. 



With this transformation, there are two unknowns at the event horizon, uj{rh) 
and (p{rh), since the field equations yield 



^ ^ ^ 1 - (1 - u^iy - + ^^)^ -m- 

where 

Uh = uirh)=uj{l) 0^ = 0(rO=0(l). (197) 

Solving the field equations (|191|) -( p!9^ ) is therefore a two-parameter shooting prob- 
lem. The procedure is to take initial 'guesses' for the unknowns Uh and (j)h and then 
integrate the differential equations out from rh using a standard ordinary differential 
equation solver, attempting to satisfy the boundary conditions for large r. The ini- 
tial starting values for Uh and (j)h are then adjusted until these boundary conditions 



are satisfied (see [B^ for further details of the algorithm used). 



For each fixed value of the Higgs mass v, there are many solutions which can 
be indexed by the number of nodes k of the potential function uj{r). Here we 
concentrate on the case k = 1. Then, for each v, there are two solutions which can 
be ascribed to one of two families of solutions: the k = 1 branch or the quasi-fc = 
branch, depending on the behaviour of the families as f —>■ 0. The quasi- /c = 
branch of solutions approaches the Schwarzschild solution uj = l,(f) = Oasv^O, 
whereas the k = 1 branch of solutions approaches the first coloured black hole of 
[ p!9[] as f — > 0, with = 0. As f increases, the two branches of solutions join up 
at t> = Vmax = 0.352. This phenomenon does not occur for A = 0.125, as found by 
Greene, Mathur and O'Neill Q. However, they conjectured that the two branches 



of solutions would converge of some value of A. We stress here that our approach is 
somewhat different from that of ref. 0, where the field equations were solved for 
fixed Higgs mass v and varying r/^, whereas we have fixed and varied v. 



For each value of the Higgs mass f , we calculated the quantities 
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{2u" + r'(l)")dr 



(198) 
(199) 



for each of the two solutions. The resulting solution curve plotted in {v, 6o, J^) space 
is shown in figure 1. The projection of this curve on to the {v,A4) plane are shown 
in figure 2. 



One issue that is important, especially when we come to consider the thermody- 
namics and entropy of the black holes, is whether or not they are extremal. An 
extremal black hole occurs when N has a double zero at the event horizon, and is 
caused physically by an inner horizon moving outwards until it coincides with the 
outermost event horizon. Mathematically, the condition for extremality is that 



m'(l) 

From the field equations ( |194| ), we have 
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(201) 



(202) 



where in the last relation we have kept an explicit Vh dependence for calculational 
convenience. There is thus no a priori reason why this quantity should not be equal 
to one half for some equilibrium solution. For the solutions on the k = 1 and quasi- 
k = branches, we can however place the following bounds on m'(l). The first term 
is decreasing for Uh positive and increasing, and hence is bonded above by its value 
for the smallest value of Uh along these branches, which is Uh = 0.632, whence 



(l-iufy < 0.360. 



(203) 



Along both these branches, 
second term. 



< 0.1 9f which gives the following bound on the 



-(plil + uohf < 2 X {Q.l9vf < 2 X 0.19^ X 0.352^ = 8.95 x 10 



-3 



(204) 



Finally, for the last term we have 

^{(pl - < ^v^ < 0.15 X 0.5 X 0.352^ = 1.15 x 10~l (205) 

Adding together all the contributions, we find that 

m'(l) < 0.5 X (0.360 + 8.95 x 10"^ + 1.15 x 10"^) = 0.185 < 0.5 (206) 
and hence all the equilibrium black holes considered here are non-extremal. 



Linear perturbation equations 

Consider small, time- dependent perturbations about the equilibrium solutions dis- 
cussed above, within the initial ansatz for the metric and matter field functions. We 
use a 5 to denote one of these small perturbation quantities, all other quantities are 



assumed to be static equilibrium functions. Following ref. [|T3[, we set 5ao = so 
that the field configurations remain purely magnetic. With this choice, the pertur- 
bation equations decouple into two independent coupled systems. The first concerns 
Sai, Suj and 6ip only. The equations take the form, with a prime denoting d/dr* 
where r* is the tortoise co-ordinate: 

- Nr'^5ai = 2N^S^ (^lu'^ + ^0^^ Sa^ + 2NS{uj6u - uj'6Cj) 

+ K'^NS{(j)'6^ - (pStp') (207) 
26i; = 2{NSuj6aiy + 2NSuj'6ai + 6Co" + NS^cpdi/j 

(208) 




= -{NSr^4>6aiy + -r^NS4>'6ai- NS^cp6Co - {r^Sipy 

+2iV^2 1^ (1-^)' + 1^2^(^2 _ ^2)^ (209) 

= +2^cD-y#^|. (210) 

This final equation is known as the Gauss constraint equation, since it represents an 
additional constraint on the field perturbations rather than an equation of motion. 
This system of coupled equations is referred to as the sphaleronic sector because it 
does not involve any perturbations of the metric functions. 



The remaining perturbation equations form the gravitational sector and concern 
the perturbations of the metric functions and also 5u! and 5(j): 

-Su = -Su" + U^Ju + U^^Sc^ (211) 
-54> = -5(1)" + U^Jtu + U^^5(l) (212) 

where the f/'s are comphcated functions of A^, S, uj and (f) and are given exphcitly 
in section 4, equation |105| . The equations governing the behaviour of 5m and 5S are 
derived from the hnearised Einstein equations and are: 

A(S,,„) . |(2iVS^,, + r=iVS^.,) (213) 

5m = 2N^5u + r^N^5^ (214) 
dr dr 



^ f I dS\ Adujd5uj d^^^^ 



From (El 31) 5m has the form 



5m = 2N^5uj + Nr^'^5ct) + ^ (216) 
dr dr S 

where f{t) is an arbitrary function of t. Compare this with the following, which 
results from integrating ( 214 ): 

5m = 2N'^5uj + Nr'^'^5<i) + g{r) (217) 
dr dr 

where g{r) is an arbitrary function of r. Comparing ( |216|) with ( |217| ), we see that 
fit) = = g{r) and 

duj r,d(J) , , 

5m = 2N—5u + Nr^-^5(p. 218 
dr dr 

We consider periodic perturbations of the form 

5cj(r,t) =(5cu(r)e^"* (219) 

and similarly for the other pertubation quantities. When substituted into the per- 
turbation equations for each of the two sectors, the equations studied in detail in 
sections 3 and 4 are derived. 



Appendix B 



Definitions and results of catastrophe theory 

Consider a family of functions 

f-.XxC^M f{x,c) = Ux) (220) 

Here X and C are both manifolds known as the state space and control space 
respectively. In other words, we have a family of functions of the variable x, the 
members of the family being indexed by c. From now on we take both X and C to 
be intervals of the real line. Then / maps out a surface z = f{x, c) in which is 
known as the Whitney surface p8| . 



The catastrophe manifold is defined as the subset of X x C at which 

^Ux) = 0, (221) 

namely it is the set of all critical points of the family of functions. In section 4, 
critical points of the functional correspond to solutions of the field equations, 
and hence the catastrophe manifold corresponds to the projection of the solution 
curve onto the (x, c) = {6o, v) plane. 



The catastrophe map x is the restriction to the catastrophe manifold of the natural 
projection 

n:XxC^C, n{x,c) = c. (222) 

This can easily be extended to a projection of the solution curve on to the (c, z) 
plane: 

Xix, c,z = f{x, c)) = {c,z = /(x, c)). (223) 

The singularity set is the set of singular points of x the catastrophe manifold, 
and the image of the singularity set in C is called the bifurcation set B. Here both 
manifolds X and C are of dimension 1, and hence x be singular whenever its 
derivative vanishes. 



The first result we require is that the singularity set is the set of points (x, c) at 
which /c(x) has a degenerate critical point, in other words, both 



dx 



fc{x) = and -^fcix) 



0. 



(224) 



This implies that the set B is the place where the number and nature of the critical 



points of the family of functions /c(x) change (see p8[ for more details of these 
results) . 



In our case, where both X and C are one-dimensional, the only possibility is that 
the bifurcation set B either is empty (in which case there is no catastrophe) or B 
contains a single point (when a fold catastrophe occurs). We observe in section 4 
that the latter situation arises. 

The catastrophe manifold is a curve C in the {x, c) plane, with the point B lying 
on it. On one side of the point B, points lying on C correspond to maxima of the 
functions fc{x) whilst on the other side of B they represent minima. In section 4 
the value of / corresponds to energy, so that minima of / will represent (relatively) 
stable objects, whilst maxima of / will represent (relatively) unstable configurations. 
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Figure Captions 



Figure 1 Solution curve for black holes in EYMH theory with one node of the 
gauge field component cu, in {v,5o,A4) parameter space: v is the Higgs v.e.v., 5q is 
the black hole parameter defined in Appendix A and Ai is the mass functional of 
the black hole. Notice that the solution curve is smooth, but this is not true for its 
projection onto the {v, M.) plane, see Figure 2 below. 



Figure 2 Projection of the solution curve of Figure 1 onto the {v, Ai) plane. The 

cusp at V = Vmax indicates the role of TVf as a fold catastrophe functional, with v 
the appropriate control parameter. The upper branch of solutions (quasi- A; = 0), 
corresponding to higher entropy, is more stable relative to the lower branch [k = 1). 
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